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We have investigated quantum phase transition employing the quantum renormalization group 
(QRG) method while in most previous literature barely entanglement (concurrence) has been demon- 
strated. However, it is now well known that entanglement is not the only signature of quantum 
correlations and a variety of computable measures have been developed to characterize quantum 
correlations in the composite systems. As an illustration, two cases are elaborated: one dimensional 
anisotropic (i) XXZ model and (ii) XY model, with various measures of quantum correlations, 
including quantum discord (QD), geometric discord (GD), measure-induced disturbance (MID), 
measure-induced nonlocality (MIN) and violation of Bell inequalities (eg. CHSH inequality). We 
have proved that all these correlation measures can effectively detect the quantum critical points 
associated with quantum phase transitions (QPT) after several iterations of the renormalization in 
both cases. Nonetheless, it is shown that some of their dynamical behaviors are not totally similar 
with entanglement and even when concurrence vanishes there still exists some kind of quantum 
correlations which is not captured by entanglement. Intriguingly, CHSH inequality can never be 
violated in the whole iteration procedure, which indicates block-block entanglement can not revealed 
by the CHSH inequality. Moreover, the nonanalytic and scaling behaviors of Bell violation have also 
been discussed in detail. As a byproduct, we verify that measure-induced disturbance is exactly 
equal to the quantum discord measured by a z for general X-structured states. 

PACS numbers: 03.65.Ud, 03.67.Mn, 03.67.-a 



I. INTRODUCTION 



The origin of research into quantum correlation, or 
more precisely, quantum entanglement can date back to 
the EPR paper [lj] in 1935, and nowadays it is no doubt 
that quantum entanglement is one of the most signifi- 
cant concepts in quantum information processing pj . It 
has already been recognized as the fundamental feature 
of quantum mechanics and utilized as a crucial resource 
for communication and computation. However, entan- 
glement should not be viewed as the unique measure 
of quantum correlations since there exist other types of 
nonclassical correlations which are not captured by en- 
tanglement. Recently many authors have proposed a va- 
riety of computable measures to characterize quantum 
correlations in the composite states: quantum discord 
(QD) ||, |j, geometric discord (GD) || ||, measure- 
induced disturbance (MID) J7j, measure- induced nonlo- 
cality (MIN) H , ameliorated MID [§ and so on. Within 
such a quantum-classical framework, a great deal of con- 
cern has been raised by quantum discord and discord-like 
correlation measures in the past few years (for review, see 
JTo[ and references therein). 

In particular, as an important application in quan- 

cntanglcment can be 



turn phase transition (QPT) |11 



exploited to determine the critical points (CP) for spin 
chains at zero temperature [0-n6|. Meanwhile, since 
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quantum discord (and other discord-like measures) is in- 
troduced as an information-theoretical tool to qualify and 
quantify quantum correlations, it is natural for us to clar- 
ify the role played by quantum discord in QPT. Several 
studies concer ning such a relationship have already ap- 
peared in Ref. |l^-|l9(j . These recent observations demon- 
strate that the CP information provided by QD is just in 
agreement with that of entanglement, and even at finite 
temperature QD still works fine Jl9[ . 

Instead of resorting to two-point spin-spin correlation 
functions, which is usually done in most previous liter- 
ature, quantum renormalization group (QRG) method 
p0| , pl| is introduced to investigate the quantum infor- 
mation properties of critical systems. Invoking such a 
method, surveys regarding Ising and Hciscnbcrg models 
have been carried out by several works [|2^-|27J and it 
has been shown that implementation of QRG method is 
valuable in detecting the nonanalytic behavior of entan- 
glement (concurrence) and the scaling behavior in the 
vicinity of CPs. Nevertheless, as mentioned above, en- 
tanglement is not sufficient to account for all the corre- 
lation contained in quantum systems, so this motivates 
us to apply other correlation witnesses to study their dy- 
namic behaviors combining with QRG method. To serve 
as a further comparison, the violation of Bell inequalities 
p8| , p9[ is also taken into consideration. 

The outline of this paper is as follows. In Sec. II and 
Sec. Ill, we investigate one dimensional anisotropic XXZ 
model and XY model respectively, using several kinds of 
correlation indicators under the method of QRG. In Sec. 
IV, we discuss the scaling behavior of these quantifiers 
when close to CPs. Sec. V is devoted to the discussion 
and conclusion. Finally, some technical points arc clari- 



2 



fied in the Appendix. 



II. CORRELATION ANALYSIS IN THE 
ANISOTROPIC XXZ MODEL 

First, we recall the QRG method and its application in 
one dimensional anisotropic XXZ model. In fact, renor- 
malization group refers to a mathematical tool that al- 
lows systematic investigation of the changes of a physical 
system as viewed at different distance scales. The key 
point to QRG scheme lies in reducing the effective degrees 
of freedom of the system through a recursive procedure 
until a mathematically tractable situation is reached. 
Following KadanofFs approach (the "block-spin" renor- 
malization group), the (one dimensional) lattice is split 
into blocks. The Hamiltonian of each block is diagonal- 
ized exactly to obtain the low-lying eigenstates (project 
operator) to construct the basis for renormalizcd Hilbcrt 
space. Finally the full Hamiltonian is projected onto the 
renormalizcd space to achieve an effective Hamiltonian 
H e ft '. Here we can summarize the QRG method as the 
following steps: 

1) Decomposing the Hamiltonian into the intrablock 
and interblock parts: H = H B + H BB , where H B is 
the block Hamiltonian, and the interblock interaction is 
denoted as H BB . 

2) Diagonalization of H B : this procedure is aiming to 
obtain the low-lying eigenstates and build the projection 
operator Pq onto the the low energy subspacc. 

3) Renormalization of H B and H BB : by virtue of per- 
turbative expansion (see Ref. [30]), the effective (renor- 
malized) Hamiltonian up to the first-order correction is 



H eff 
T eff 



H, 



T eff 



where H, 



eff _ 



PnH B Pn and 



HI 1 3 = P H BB P . 

4) Iteration: repeat 1) 3) to arrive at the final man- 
ageable situation. For more details, we refer the readers 
to Ref. [30-32]. 

Kargarian ct al. introduced the notion of "renormal- 
ization of concurrence" ^2j, and they found that this 
notion truly captures the nonanalytic behavior of the 
derivative of entanglement (concurrence) close to the crit- 
ical point. As a warmup, we first review the renormaliza- 
tion of entanglement in the one-dimensional anisotropic 
XXZ model (23). Thc Hamiltonian of spin 1/2 XXZ 



model on a periodic chain of N sites is 



N 



H(J, A) = -J2 Kof+i + TjVlVi + Aof of +1 ) , (1) 

i 

where J is exchange constant, A is the anisotropy pa- 
rameter, and J, A > 0. erf (a = x,y,z) are standard 
Pauli matrices at site i. This model is known to be 
exactly solvable by Bethe Ansatz and critical (gapless) 
while < A < 1. The Ising regime is A > 1 and a 
maximum of concurrence can be reached between two 
nearest-neighbouring spins at the transition point A = 1 

MM. 



To construct a renormalizcd form for thc Hamiltonian 
(|l|), we shall choose a decomposition of three-site blocks. 
Note that this is requisite in the sense that it is a guar- 
antee of self-similarity after each iterative step. Ref. ||^| 
gives the degenerate ground states of the block Hamilto- 
nian as follows 



|0o 



1 



=(im) + gint) + utt», (2) 



\4>' ) = -^=(i m> + q\ m) + 1 ut)), (3) 

where | f) , I I) ar e the eigenstates of er 2 and 



<Z = -^(A+ VA 2 + 8), 



(4) 



Then the effective Hamiltonian of thc renormalizcd chain 
can be cast into thc form 



, N/3 



where the iterative relationship is 



16) 



The most important information given in the QRG 
method are its fixed points. By solving equation A' = A, 
we obtain the trivial fixed point A = and also the non- 
trivial fixed point A = 1. It is worth noticing that, as was 
stated previously, XXZ model is critical for all < A < 1 
while QRG method only indicates the single points. In- 
deed, if appropriate boundary terms are implemented in 
the QRG method, then it predicts correctly a line of crit- 
ical models in the range < A < 1 |j35| . 

In order to calculate quantum discord and other cor- 
relation quantities, we consider one of the degenerate 
ground states. Correspondingly, the density matrix is 
defined as 



P123 = |<M(<A)| 



(7) 



with \(j>o) referring to Eq. (g) (choosing \<f)' ) will yield the 
same results). Since we are focusing on pairwise corre- 
lation functions, without loss of generality, we trace over 
site 2 to obtain the reduced density matrix between site 
1 and 3 



Pl3 



( q 2 

110 

110 

V o o o o 



(8) 



It is straightforward to compute the concurrence [^6| of 

Pl3 



c 



13 



(9) 
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It is shown that the concurrence between two blocks ex- 
hibits an explicit signature of the quantum phase transi- 
tions at A = 1. Meanwhile, we notice that some other en- 
tanglement measures in the literature, such as von Neu- 
mann entropy or the averaged bipartite entanglement, 
are shown to be good indicators of the quantum phase 
transition. However, to our best knowledge, there is no 
study of analyzing other correlation witnesses beyond en- 
tanglement in QRG framework. In the below section, we 
will analytically derive these quantities in detail to see 
whether they can be proved really helpful in predicting 
critical phenomenon. 



Quantum discord and measure-induced 
disturbance 



a two-qubit X-shaped state. A general X-state have the 
visual appearance 



( pu p 14 

P22 P23 

P*23 P33 

\ p* u /9 4 4 



(16) 



which has seven real parameters. However, up to local 
unitary equivalence, we can assume pi4 and P23 are also 
real and in fact there are only five independent parame- 
ters (note that QD is invariant under local unitary trans- 
formations). Alternatively, if we represent the X-state 
in Bloch decomposition, then the five characterizing pa- 
rameters can be expressed as 



Quantum discord is introduced in Rcf. [|3j aiming to 
characterize all the nonclassical correlations present in a 
bipartite state. It originates from the inequivalence of 
two expressions of mutual information in the quantum 
realm. Consider a composite bipartite system p AB , the 
quantum mutual information is defined as 



x 

II 
h 
fa 
f 3 



Tr(a A p x ) = pu + P22 
Tr(a B p^ - 



P11 - P22 



P33 
■ P33 



P44, 
P44, 



Tr(a A a^) = 2 Pli + 2p 23 , 
Tr(a A a B p x ) = -2p 14 + 2p 23 , 

Tr(a A a B p x ) = p X \ - P22 - P33 + Pa, 



(17) 



T{p A »):=S{p A ) + S{p»)-S{p A ») 



(10) 



where S(p) = —Tr(p\og 2 p) is the von Neumann entropy, 
and p A ( B > = Ttb(a)(p AB ) denote the reduced density 
operator of subsystem A(B). On the other hand, if a 
complete set of von Neumann measurements {n^} (or 
more generally, POVMs) performed on subsystem A, an 
alternative version of quantum mutual information con- 
ditioned on this measurement yields 

l(p\{U A }): = S(p B )-S(p\{U A }), (11) 
= S(p B ) -J2PkS(p B ), (12) 

k 

with Pk = Tr(U A p AB ) and p B = Tr A (U A p AB )/ Pk . To 
eliminate the dependence on specific measurement, one 
takes the optimization procedure to obtain 



Except for some numerical evaluations for a restricted 
subset of two-qubit X states |37], [$8|, Ref. Q has pre- 
sented a algorithm to calculate QD for general two-qubit 
X states: the optimal measurement is in a universal finite 
set {a x , a y , a z } (see also Ref. p0|, j4l| ). Nonetheless, a 



counterexample is given by Rcf. §42J to disprove the al- 
gorithm and this fact elucidates the state-dependence in 
the optimization for general X states. Furthermore, re- 
cent progress toward this problem has been made by Ref. 
p0| , |4l[| , which identifies a large class of X states whose 
QD can be derived analytically from the measurement 
strategy described above. 

Keeping these technical preparations in mind, let's 
turn to density matrix (B|), the spectrum of which is 

KP13) = {0, 0, 2+^2-}- First we note that p 22 
and then 



P33 



J{p) :=m^l(p\{Hf}), 



(13) 



which has been suggested by Henderson and Vedral ||] 
as a measure to quantify the purely classical part of cor- 
relations. The discrepancy between the original quan- 
tum mutual information X and the measurement- induced 
quantum mutual information J is defined as the so called 
quantum discord 



V A (p):=l(p)-J(p), 

= S(p A ) - S(p AB 



) + min S2p k S(p 
{n^} V 



(14) 
(15) 



Now we are going to deal with the situation that we 
come across. The density matrix defined in Eq. (H) is 



S(p A ) = S(p B ) 



q- 



(18) 



which means T>a{p) = T^b{p) U (here A and B denotes 
site 1 and 3 respectively). According to the Theorem in 
Ref. [^0|, it is easy to check that the optimal observable 
for state (||) is not a z but a x . For the reason that will 
be clear later, we also gives discord measured by a z (for 
brevity we put the general formulas in the Appendix) 
which is actually equal to the concurrence 



A 



(19) 
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When the optimal measurement {-^(Ha A )} is performed 
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on A, wc directly acquire the QD of state 
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(20) 



where /(z) := -±^log 2 - 



l^]o ga i^. Nu- 
merical evaluation shows 2?^ = 2?^ is strictly less than 
, as we expect. 
Based on the definition introduced by Luo |7j , measure- 
induced disturbance (MID) is defined as the difference of 
quantum mutual information before and after measure- 
ment 
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FIG. 1: (Color online) The evolution of the quantum discord 
versus A in terms of QRG iteration steps in the XXZ model. 



MID(p AB ) = I(pab) - HTI(pab)), 



(21) 



The measurement II 



{n^(g)II^} is induced by the spec- 



tral decompositions of the reduced states, p A = Y^pip A 
, which leaves the marginal infor- 
By taking ttf = and Uf = 

(which is unique in our case), we have XI(pab) = 
diag{q 2 , 1, 1, 0} and consequently 



and p B = T,iPip 
mation invariant. 



1 

2+q 1 



MID = S(U(p AB )) - S(pab) = 



(22) 



So MID of state (||) coincides with T> a £ and concurrence. 
However, this is not a coincidence. In the Appendix we 
will prove that MID is exactly equal to the QD measured 
by a z for general X-structured states. 

In Figure [j], we illustrate the evolution of QD versus A 
for different QRG steps. Notice that the iterative rela- 
tionship we adopt here and later in the calculation is Eq. 
(4) and (6). The plots of QD cross each other at the criti- 
cal point A = 1. In comparison with concurrence demon- 
strated in Ref. [^3| , QD also develops two saturated val- 
ues, which are associated with the two different phases, 
spin-liquid and Neel phases. Note that after enough it- 
eration steps for < A < 1, D A m 0.412154 < C = 0.5 
while for A > 1, Da —> 0. In addition, since MID is 
exactly equal to concurrence, MID can obviously exhibit 
a QPT at A = 1. 



Geometric discord and measure-induced 
nonlocality 



For two-qubit case, a closed form of expression for geo- 
metric discord can be achieved 



\T\\ 2 



(24) 



where Xi = Tr(a A p) are components of the local Bloch 
vector for subsystem A, Tij — Tr{a A a B p) arc compo- 
nents of the correlation matrix, and x := (x\,X2,Xs) , 
T := (Tij), X m ax is the largest eigenvalue of the matrix 
K = xx^+TT* (here the superscript t denotes transpose). 
It is worth emphasizing that, Luo and Fu presented an 
equivalent but simplified version of the geometric discord 



mm||p-n^(p)|| 2 , 

\1 A 



(25) 



where the minimum is over all von Neumann measure- 
ments II = {n^} on subsystem A. Intuitively, in some 
sense dual to GD, another measure quantifying nonlo- 
cal effect caused by locally invariant measurements was 
introduced Luo and Fu [|j 



MINa(p) =max||p-n A (p)||' 



(26) 



with an extra constraint that von Neumann measure- 
ments H A = {n^} do not disturb p A locally, which means 

/ = E t nyn^. 

For X state and its characterizing parameters de- 
fined in Eqs. ([l7|), we have x = (0,0, a;)' and T = 
diag{ti,t2,t^} and GD of X state reads 



In this subsection we calculate GD and MIN. Recently, 
Dakic et al. introduced the following geometric measure 
of quantum discord based on the Hilbcrt-Schmidt norm 



1 



V^^ti + ti+ti + x 2 
- max{t 2 ,^,^ 



x 2 }), (27) 



V G A (p) :=min||p- X | 



(23) 



where denotes the set of zero-discord states and ||p — 
x|| 2 = Tr(p~x) 2 is the square of Hilbert-Schmidt norm. 



Besides, using Theorem 3 in Ref. 
tained for two-qubit X states 



MIN a (p x ) 



MIN can be ob- 



if x ^ 



(tl + q+q- \ mm ), if x = 



(28) 



5 



with X miT 
the state 



= min{f^, t\, if.}. Applying these formulas to 
]), we obtain 



1 



V<l = MIN A = t(tl+tl) = 



by noting that t 



,3 - 



>-2+a 2 / +t -2+o 2 / — l >2+S 2 "'' 



t 



t\ and a; 



2+g 2 > 1 V2+q 2 ^ — V 2+q 2 ■ 

since |g| > \/2 for A > 0. The variation 
of GD and MIN versus A has been plotted in Figure |[ It 
is no surprise that they can indicate the precise location 
of the critical point A = 1 since T>^ = MIN A = \C 2 . 



where Ui, i = 1,2,3 are the eigenvalues of U — T f T. 

As to X states, the matrix T is diagonal and T*T = 
diag{t\ , if. Therefore the maximal violation of the 
CHSH inequality for X states can be simplified to 



(33) 



with X m in = min {t%, i|, tg}. Thus, the (maximal) Bell 
violation of state (0) is obtained 
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FIG. 2: (Color online) Fhe evolution of the geometric discord 
and measure-induced nonlocality versus A in terms of QRG 
iteration steps in the XXZ model. 
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C. Bell violation 

Quantum nonlocality, as revealed by the violation of 
Bell-type inequalities, refers to many-system measure- 
ment correlations that cannot be simulated by any local 
hidden variable theory. In particular, for two-qubit pure 
states, the presence of entanglement guarantees violation 
of a Bell inequality (Gisin's Theorem) p3[ . However, 
for mixed stares the situation becomes more complicated 
p4j . Here we restrict ourself to the CHSH inequality. 
The Bell operator corresponding to CHSH inequality can 
be formulated in the following form 

Bchsh = a ■ a ® (b + 6') • a + a' ■ a (8 (6 - b') ■ <rj(30) 

where a, a', b, b' are the unit vectors in R 3 , and cr = 
(a x , <J y , cr z ). Then the well-known CHSH inequality is 
expressed as 

B = \(Bchsh) p \ = \Tr(pB C HSH)\ < 2. (31) 

According to the Horodecki criterion pq [, the maximum 
violation of CHSH inequality is given by 

B^hsh = max Tr(pBcHSH), 

a, a ,o,o 

= 2 /max(ui + %■)> (32) 

V »<? 



FIG. 3: (Color online) The evolution of the violation of CHSH 
inequality versus A in terms of QRG iteration steps in the 
XXZ model. 

The value of maximum Bell violation under QRG it- 
erations is displayed in Figure [| Interestingly, we can 
observe that the block-block correlations never violate 
the CHSH inequality, but still evidently exhibit a QPT. 
At the critical point A = 1, the (maximum) Bell viola- 
tion is a fixed nonzero constant B c = 2^2/3 w 0.942809, 
irrespective of the iterative steps. What's more, beyond 
the critical point, Bell violation also develops two satu- 
rated values: one is B = y/2 for < A < 1 and one 
is B = 2 for A > 1, which is in sharp contrast to the 
behavior of concurrence. 



III. CORRELATION ANALYSIS IN THE 
ANISOTROPIC XY MODEL 

In this section we embark on studying the relationship 
between the QPT and quantum correlation witnesses in 
the spin-1/2 XY model with the QRG method. The 
Hamiltonian of the XY model on a periodic chain with 
N sites reads 

T N 

H ^ 7) = 4 E K 1 + + (! - 7KVf +1 ] , (35) 



6 



where J is exchange coupling constant, 7 is the 
anisotropy parameter. The XY model reduces to the 
XX model for 7 = or the Ising model for 7 = 1. In 
the parameter range < 7 < 1, it falls into the Ising 
universality class. To implement the QRG method, we 
still choose three sites as a block. In Ref. |^6| Ma et al. 
obtained two doubly degenerate eigenvalues of the block 
Hamiltonian as below 

1*0) =— =L=(-v/it^i m> + V2| m) 

-Vl+flltt) + V2 7 \iU)), (36) 

i$' ) =_L=(-V2 7 | nt> + m) 

2^1 + 7 2 



V2\ m) + v/r+T^i wt», 



(37) 



After projection onto the renormalized subspace, the ef- 
fective Hamiltonian can be written as 



J' 



N 



HCff = T E K 1 + + (1 - , (38) 



with the iterative relationship 

3 7 2 + 1 , 7 3 + 3 7 



J' = J 



2(1 + 7 2 )' 7 3 7 2 + l' 



(39) 



Naturally we are most concerned with the CP informa- 
tion. The stable and unstable fixed points can be gotten 
by solving 7' = 7. The stable fixed points locate at 
7 = ±1, and the unstable fixed point is 7 = which 
separates the spin-fluid phase (7 = 0) from the the Neel 
phase (0 < | 7 | < 1). 

Similarly, we consider one of the degeneracy ground 
states to construct the pure-state density matrix 



0123 = |$o)($o|, 



(40) 



The result of choosing |$q) will be the same. By tracing 
out site 2, we arrive at the reduced density matrix 



2 2 7 

7 2 + l 7 2 + l 

013 4( 7 2 + l) 7 2 + l 7 2 + l 

2 7 2 7 2 



(41) 



The concurrence between the sites 1 and 3 is given as 
1 l7l 



2 I + 7 2 



(42) 



The reason is twofold: bringing in such a variable is not 
only convenient for us to compare the results with that of 
p6[ , but also useful in the derivation process. According 
to Ref. [|o[ [tH , h can be verified that a z is not the 
optimal observable and the choose of optimal observable 
depends on the value of 7, or more accurately, the sign of 
7: the the optimal observable to measure is <r x if 7 > 
(which is equivalent to t\ > t<i or |g| > 1), and a y if 7 < 
(|g| < 1). In spite of this fact, we still provide the discord 
measured by a z here 
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2( 7 2 + l) 82 2( 7 2 + l) 
log: 
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(44) 



We again remark that 022 = £33 for the density matrix 
@ and thus 



S( P A ) = S(p B ) = 



° (45) 



4( 7 2 + 1) V 3 7 2 + 



So we do not need to specify on which subsystem the 
measurement is performed. The spectrum of ([4l| ) is 
Mq) = {1/2: 1/2,0,0}, irrespective of the value of 7. Af- 
ter some lengthy but standard algebra, one finally gets 
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(46) 



with y, t\, t-x and function /(•) defined the same as above. 

From the proof in the Appendix, we already know that 
MID is exactly equal to V a *. The QD and MID be- 
tween site 1 and 3 have been plotted in Figure || (recall 
that for this model the iterative relationship is just Eq. 
(39)). The dynamic behavior of QD in each iteration step 
is analogous to that of concurrence, but not totally the 
same. At the critical point g = 1, QD reaches a nonzero 
constant T> m 0.412154, which once again indicates that 
spin-fluid phase contains quantum correlations as already 
shown in XXZ model. In contrast to QD, MID also shows 
the nonanalytic property, however, for < g < 1 and 
g > 1 MID does not fall to zero but gets to another 
nonzero constant MID = 1. 



A. Quantum discord and measure-induced 
disturbance 

Before calculating QD and other correlation quantities, 
we regard all these measures as a function of g, where 



B. Geometric discord, measure- induced nonlocality 
and Bell violation 

In section II we have obtained the analytic formulas 
for GD, MIN and Bell violation. Employing Eqs. @, 
(||) and (||), the GD, MIN and Bell violation for the 
state (El|) are listed as follows (in the derivation note 
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FIG. 4: (Color online) The evolution of quantum discord and 
measure-induced disturbance versus A in terms of QRG iter- 
ation steps in the XY model. 



FIG. 5: (Color online) The evolution of the violation of CHSH 
inequality versus A in terms of QRG iteration steps in the XY 
model. 



relations between the two renormalized sites represents 
the correlations between two parts of the system each 
containing N/3 sites effectively. In this sense, we can 
refer to these quantities considered in this work as block- 
block correlations. It has been demonstrated that the 
first derivative of all these correlation measures shows a 
nonanalytic behavior in the vicinity of the critical point . 
Furthermore, the scaling property of entanglement has 
also been observed in Rcf. [£3|, pq , which is related to 
the divergence of the correlation length as the critical 
point is approached. 



that 7 > and 7 < correspond to |g| > 1 and |g| < 1 
respectively) 



MIN 

Bm 



- - 

4 \ 2 

7 



7 



= 1 c 

I + 7 2 ; 4 ' 



4 ' 6 7 2 + l 



8( 7 2 + l) 2 



12 7 2 



7 2 + l 



Wmin. 



(47) 
(48) 
(49) 



Here the Bell violation versus g changing for different 
iterations is depicted in Figure [| It is clearly seen that 
in XY model the block-block correlation is still not strong 
enough to violation the CHSH inequality although Bell 
violation displays the nonanalytic behavior at g — 1 (7 — 
0). 




FIG. 6: (Color online) The logarithm of the absolute value of 
minimum, \n(\dB / dg\ m ) , versus the logarithm of chain size, 
ln(iV), which is linear and displays a scaling behavior (B is 
the Bell violation in the XY model). 



IV. NONANALYTIC AND SCALING 
BEHAVIOR 

So far, we employ the QRG method to investigate the 
block-block correlations of one-dimensional XXZ and XY 
spin models. As we have described in the QRG approach, 
the size of a large system (N = 3™ +1 ) can be effectively 
rescalcd to three sites with the renormalized couplings 
of the nth RG iteration. Therefore, the quantum cor- 



Aiming to compare with the previous results concern- 
ing entanglement, wc explicitly show the nonanalytic 
phenomenon and scaling behavior of other correlation 
witnesses. Here, we focus on the dynamic property of 
nonlocality, that is, the maximal violation of Bell-CHSH 
inequality, since our results display that quantum dis- 
cord behaves more similarly to entanglement (see Figure 
and |). First, we have analyzed the scaling behav- 
ior of y = \dB/dg\ glll versus the size of the system N 
in the XY model where g m is the position of the mini- 



8 



mum of dB/dg. We have plotted ln(y) versus ln(iV) in 
Figure ^ which shows a linear behavior. Numerical cal- 
culation tells us that the exponent for this behavior is 
\dB/dg\ gxn ~ N°". As a companion, our analysis also 
reveals that the position of the minimum g m of dB/dg 
gradually tends to the critical point g c = 1, as shown 
in Figure ^ (the numerical relation is g m = g c — N -1 ' 00 ). 
These results convince us that Bell violation truly signify 
the criticality of the spin system. 



-2 




2 4 6 8 10 12 



ln(N) 



rithm raised in Ref . |40j, |4l| , where the optimal measure- 
ments to achieve QD can be exactly determined. Besides, 
we are convinced that the whole analysis in this paper 
can be extended to many other spin models, since the re- 
duced density matrices are usually highly symmetric and 
can be casted into X-shaped states |24| , |25| , p7| , |3l[ [32| . 
Very recently, it has been reported that Bell inequality 
is able to signal QPT and it can never be violated in the 
corresponding spin models jl6[ |47j. In these works, the 
nearest-neighboring spin-spin correlation functions are 
invoked to compute the Bell violation, which are usually 
complex and lengthy. However, we resort to the QRG 
framework and also illustrate no violation can be discov- 
ered in each iteration step, which implies some intrinsic 
feature of long-scale corrections. The connection between 
these observations will be attractive and may need fur- 
ther investigation. Finally, we would like to mention that 
it might be interesting to apply the same approach to 
high-dimensional systems, where a straightforward nu- 
merical analysis could be performed for some measures 
of quantum correlations. 



FIG. 7: (Color online) The scaling behavior of g max in terms 
of system size N where gmax is the position of the minimum 
derivative of Bell violation in the XY model. 



V. DISCUSSION AND CONCLUSION 

In this paper we have investigated the performance 
of various correlation measures in quantum phase tran- 
sition, exploiting the quantum renormalization group 
method. In most previous literature along this line, 
only entanglement (concurrence) has been utilized as an 
information-theoretic tool to evaluate the critical prop- 
erties of the spin systems. However, it is now well known 
that entanglement can not account for all the aspects of 
quantum correlations, which in turn motivates us to clar- 
ify whether other correlations witnesses (including Bell- 
CHSH violation) are useful in such a circumstance. In- 
deed, there are several points that deserves our attention: 
(i) The quantum discord and other discord-like measures 
turn out to be as good as entanglement to detect the 
quantum phase transition in the anisotropic XXZ and XY 
models. Nevertheless, it is apparent that the dynamic 
processes of these quantities arc not totally similar with 
entanglement and even when concurrence vanishes there 
still exists some kind of quantum correlations which is not 
captured by entanglement, (ii) Interestingly, our result 
shows that CHSH inequality can never be violated in the 
entire iteration procedure, which indicates block-block 
entanglement can not revealed by the CHSH inequality. 
Moreover, the nonanalytic and scaling behaviors of Bell 
violation have been justified by numerical calculations. 

On the other hand, the two cases handled in this work 
can be regard as perfect examples to apply the QD algo- 
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Appendix A: Analytic proof of MID — V* for 
general X states 

Here, wc show that measure-induced disturbance is ex- 
actly equal to the quantum discord measured by a z for 
general X-structured states. The expressions of QD and 
MID arc formulated as follows 

V A {p) = S{ P A ) ~ S{p) + min y> fc S(p£), (Al) 

MID(p) = I( P ) - I(U(p)) = S(Tl(p)) - S(p), (A2) 

Note that the two-side measurements n = {Ilf ® ^f} 
employing in MID depend on the spectral decompositions 
of the reduced states. Since we only consider the discord 
measured by a z , that is, {§(/ ± a z )} = {|0)(0|, |1)(1|}, 
the conditional states for general X states (O) can be 
obtained 

P ° P11+P22 V P22 J' 

P? = ~ J — ( P n ° ), (A3) 

P33 + P44 V P** J 
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with p = pn + p 2 2 and p x 
states are all diagonal states 



P33 + PiA- The reduced 



P 



P11 + P22 

P33 + P44 

P11+P33 \ 

p 22 + P44 J ' 1 J 

Therefore, we can take H A = Hf = = 

0,1). To prove MID = V 7 *, all we need is to verify 
S(P A ) + Y, k PkS(Pk) = S(n(p)). In fact, it is easy to 
find 



s(p A ) + "£p k s(p 



-^2 Pit l0g 2 (/9«) 

ii 

S(n(p)) (A5) 



with m = 11,22,33,44. Moreover, when the measure- 
ment a z is performed on subsystem B, the situation is 
the same. To sum up, we arrive at the relationship 
MID for X states. 



n 



B 



Appendix B: Analytic formula of © <Tx (2? CT «) for 
general X states 



with p k = Tr(n A <g) I B pH A ®I B ) = \,k = ± and y, h 
defined in Eqs. ([l?]). In addition, p± have exactly the 
same spectrum 



(B2) 



Therefore 



where f(z) 



Y,PkS(p B ) = S(p B ) = S( P B ), 

k 

= fiv 2 + t\). 



^62 I 1°S2 



(B3) 



If we 



choose {11^ = |(/ ± (Tj,)} as the local measurement on 
A, an analogous expression can be achieved 



J2p k S(p B )^f(y 2 + t 2 2 ), 

According to Ref . ffl, |4l[] , as long as 



(B4) 



By definition, we need to evaluate the conditional state 
p B and the corresponding probability pk, since S(p A ) and 
S(p AB ) are easy to compute. Let {11^ = i(I ± cr x )} 
(fc = ±) be the local measurement for subsystem A, then 



p B = —Tr A (n A ®I B pn A t 
P± 

P11 + P33 ±{pu 
±(pu + P23) P22 H 

1 / 1 + y ±ii 

2 V ±*i l-y 



- P23) 

P44 



(Bl) 



\y/pllPU - \7P22P33l < \Pli\ + \p2a\, (B5) 



holds, the optimal observable is a x if ti > t,2 and (Ty 
otherwise. Consequently, 



■D A (p) = S(p A ) - S(p AB ) + f(y 2 + max{t? ,}), (B6) 



Note that S(p A ) and S(p AB ) can also be represented by 
parameters defined in Eqs. (|l7|) (sec Eq. (8) in Ref. pl[). 
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